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We study the magnetic susceptibility of one-dimensional S = 1/2 antiferromagnets containing 
non-magnetic impurities which cut the chain into finite segments. For the susceptibility of long 
anisotropic Heisenberg chain-segments with open boundaries we derive a parameter-free result at 
low temperatures using field theory methods and the Bethe Ansatz. The analytical result is verified 
by comparing with Quantum-Monte-Carlo calculations. We then show that the partitioning of the 
chain into finite segments can explain the Curie-like contribution observed in recent experiments on 
Sr2Cui_ a; Pd a; 03+5. Possible additional paramagnetic impurities seem to play only a minor role. 
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Measurements of the magnetic susceptiblity, Xi 011 an ~ 
tiferromagnetics at low temperatures, T, often reveal 
a (Curie) term \ l/T. This is usually associated 
with paramagnetic impurities (free spins). For one- 
dimensional S=l/2 systems, however, a second mecha- 
nism exists which leads to a contribution with somewhat 
similar T-dependence. Non-magnetic impurities can cut 
the chain into finite segments with essentially free bound- 
aries. At T/J < 1/L (J being the exchange constant) 
a segment with odd length L will lock into its doublet 
ground state yielding also a 1 /T contribution [1-4] . How- 
ever, the behavior is considerably more complicated at 
higher T. As we show, fitting data to a Curie form for 
a limited range of T can lead to an underestimate of the 
impurity concentration by as much as a factor of ten! In 
general, the different T-dependence could make it pos- 
sible to distinguish this type of impurity from paramag- 
netic ones providing useful structural information. 

These results are relevant for materials like S^CuOa+a 
which is known to be an almost ideal realization of the 
spin-1/2 Heisenberg chain with a nearest-neighbor cou- 
pling constant J ~ 2200 K and a low Neel temperature 
Tjv ~ 5 K~ 0.002 J [5, 6]. Measurements of the magnetic 
susceptibility have revealed a Curie contribution which 
could be dramatically reduced by annealing. The Curie 
term therefore is believed to be mainly caused by ex- 
cess oxygen [5, 6]. Recently, susceptibility measurements 
have also been reported on Sr 2 Cui_ a; Pd a: 03 + 5 with Pd 
serving as a non-magnetic impurity [7], 

In this letter we show that these susceptibility data 
can be explained by taking the segmentation of the chain 
due to the non-magnetic impurities into account properly 
without including any additional paramagnetic impuri- 
ties. To this end we derive a parameter-free analytical 
result for the susceptibility of an anisotropic Heisenberg 
chain with finite length and open boundary conditions 
(OBCs) at low temperatures using field-theory methods. 
This extends and generalizes previous results in the scal- 



ing limit [2, 3] and for the boundary susceptibility in the 
thermodynamic limit [8-11] . We will verify our analytical 
result using Quantum-Monte-Carlo (QMC) calculations. 

We consider an anisotropic spin-1/2 Heisenberg chain 
consisting of segments with length L and OBCs 
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Here, A parameterizes the exchange anisotropy. We as- 
sume that T Tjv so that interchain couplings can 
be safely ignored. As argued at the end of this letter, 
longer-range interactions which could bridge between the 
chain segments can be neglected as well in the parameter 
regime we are interested in. For a concentration of chain 
breaks p the average chain length is given by L = l/p — 1 
and the averaged susceptibility by [12] 



= P 2 J2 L (1-P) L X(L) 



(2) 



We now calculate the susceptibility x(fi) = 
((Y,iS?) 2 )/(LT) in the low T large L limit using field- 
theory methods. We start with the Hamiltonian (1) in 
the scaling limit, i.e., ignoring irrelevant operators. The 
Hamiltonian is then equivalent to a free boson model [13] 



H 



dx [n 2 + {d x <t>) 2 



(3) 



Here </> is a bosonic field obeying the standard commu- 
tation rule [(f>(x),H(x')] = iS(x — x') with n = v~ 1 dt(f>. 
The velocity v is a known function of the anisotropy A. 
Using the mode expansion for OBCs [2] 
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where R is the compactification radius of the bosonic 
field, the Hamiltonian (3) can also be expressed as 



H = J7T S * ~ hS * + IT E n ( fl « a " + V2) • (5) 



Here a„ is a bosonic annihilation operator, 6" 2 an integer 
(half-integer) for L even (odd) and h the magnetic field. 
K - the so called Luttinger parameter - is related to the 
compactification radius by K = 1/(2ttR 2 ) and is a known 
function of anisotropy A as well. The susceptibility in the 
scaling limit is then given by 
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For LT jv — ► and L even Xs ~ 7^ exp [— ^2t ] whereas 
for L odd %s ~ (4LT) -1 . For LT/t) — > 00 the thermo- 
dynamic limit result Xs — K/ (2ttv) is recovered. Correc- 
tions to scaling occur due to irrelevant bulk and bound- 
ary operators. The leading bulk irrelevant operator for 
< A < 1 is due to Umklapp scattering yielding the 
following correction to (3) 



5H = Ai / dxcos(2<p/R) 
Jo 



(7) 



with Ai being the Umklapp scattering amplitude. To 
obtain \ to first order in the operator (7) for L and 
T finite the expectation value (exp(±2i</>/i?)) has to 
be calculated. Using the mode expansion (4) this 
expectation value splits into an S z (zero mode) and 
an oscillator part. Upon using the cumulant theo- 
rem for bosonic modes we obtain (exp(±2i^>/i?)) = 
(cxp(±2i(/>/i?))s 2 exp(— 2{(j>(j>) osc. /R 2 ) with 
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Introducing a cutoff a for the zero temperature part in 
(8) and using Ei^i = — ln(l — z) for \z\ < 1 we 
obtain the following correction to (6) 
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Here r)(x) is the Dedekind eta- function, 8i(u, q) the ellip- 
tic theta function of the first kind [14] and Ai = a 2K X\. 



For K < 3/2 the integral in (9) is convergent, the cutoff 
can be dropped, and the last line becomes equal to one. 
An additional correction arises due to the presence of an 
irrelevant boundary operator. In essence, this operator 
changes the length of the chain L to some effective length 
L + a [2] . Replacing L by L + a in the exponentials in 
(6) and expanding to lowest order in a yields 



5X2 - -TTva gi (exp[-TTv/(KLT)})/(KT 2 L 3 ) (11) 
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For K > 3/2 we isolate the cutoff dependence by sub- 
tracting the Taylor expand of the integrand in (9) to first 
non-vanishing order in y, and then take a — > 0, yielding 
a convergent integral. The cutoff dependent part has ex- 
actly the same form as (11). We can therefore think of 
a in (11) as a parameter incorporating both contribu- 
tions. In the thermodynamic limit 171 (e^ 7TV ^ KLT ^ — > 
K 2 T 2 L 2 /(2tt 2 v 2 ) and 5 X 2 -» Ka/{2nvL). In this limit 
we can compare the field theory result with a recent 
calculation of the boundary susceptibility based on the 
Bethe ansatz [9, 10] leading to 

a = 2- x < 2 sin [nK/(4K - 4)] / cos [n/(4K - 4)] . (13) 

Because the amplitude Ai is known exactly as well 
[10, 15], the obtained result is parameter-free. In Fig. 1 
a comparison between this field theory result and QMC 
data for A = 3/4 is shown and in the left inset the correc- 
tions to the scaling limit (6) due to the leading irrelevant 
bulk and boundary operator are visualized. More gen- 
erally, we always find excellent agreement if T / J < 0.1 
and L > 10. The range of validity can be extended in 
principle by including terms in higher order perturbation 
theory in Ai . However, if T/v 3> 1/L it is justified to sep- 
arate x m to a bulk, x bulk , and a boundary contribution, 
Xb, and to perform the thermodynamic limit for each of 
them separately. A combination of this result [10] (blue 
dashed line in Fig. 1) with the result derived in this let- 
ter yields an excellent parameter-free description of the 
susceptibility for anisotropic chains of finite length up to 
T-J/2. 

For K = 1 (A = 1) perturbation theory alone is not 
sufficient because Umklapp scattering becomes marginal. 
In this case the coupling constant Ai has to be replaced by 
a renormalization group (RG) improved running coupling 
constant [15, 16]. Ai then becomes a function of the 
two length scales L and v/T. For small enough energies 
the smaller scale will always dominate and the running 
coupling constant becomes Xi(L, v/T) — g/A with [15] 
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FIG. 1: (Color online) QMC data for A = 3/4 and L — 99 
(red squares), L = 100 (black dots) compared to the field 
theory result derived here (green solid lines) . The blue dashed 
line denotes x = Xbuik + Xb/L (see Eqs. (148, 149) in [10]). 
Left inset: Corrections to scaling for even length chains and 
A = 3/4. QMC data are represented by symbols, the field 
theory by the green lines. Right inset: Same as main figure 
for A = 1 with a = 5.8 and Xbuik , Xb as obtained in [8, 15]. 

where 7 is the Euler constant. S\i in the isotropic case is 
therefore given by (9) with the last line dropped and Ai 
replaced by Xi(L,v/T). Furthermore, K in (6) should be 
replaced by l+g(L,T)/2 so that KL in the exponentials 
gets replaced hy L(l+ g/2 + a/L). For this case a is not 
known and we will use it as a fitting parameter. This way 
we obtain excellent agreement if T/J < 0.1 and L > 10 
(see right inset of Fig. 1). At LT <C v the corrections to 
the scaling formula (6) in the isotropic case due to the 
marginal interaction can also be understood as follows. 
In this limit the most important correction to (6) arises 
from the correction to the excitation energy of the lowest 
excited states with S z = ±1, E -> (wv/L)[l - g(L)], [17] 
corresponding to the replacement 1/K — > l — g{L) in (6). 
Here g{L) is given by (14) but with min[L, v/T] simply 
replaced by L. 

These analytical results, supplemented by numerical 
data for L < 10, can be used to calculate the averaged 
susceptibility (2) at low temperatures. In addition, we 
have calculated x(£) numerically for all L G [1, 100] and 
certain additional L up to L = 1000 at various tempera- 
tures. An interpolation and extrapolation was then used 
to obtain x(L) for all lengths up to L = 5000 at a given 
temperature [19]. This allows us to reliably calculate \p 
(2) for impurity concentrations p > 0.2% (L < 500). 

Experimental data for the susceptibility of the quasi 
one-dimensional compounds Si^Cui-^Pd^Os+a have 
been analyzed in [5-7] by decomposing it into x = 
Xo + Xc + Xint. • Here xo represents a constant part due to 
core diamagnetism and Van Vleck paramagnetism, xc a 
Curie term and Xinf. the susceptibility (per unit length) 
of the infinite 5=1/2 Heisenberg chain. However, our 
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FIG. 2: (Color online) T Xp for p = 0.01,0.03,0.05,0.07,0.1 
from bottom to top (black dots with dashed lines being a 
guide to the eye). The blue solid lines denote the value (p/4) • 
(1— p)/(2 — p) reached in the limit T <C pv. Inset: Same data 
for Tx P with subsequent curves shifted by 0.01 compared to 
TXbuik + p/[121n(2.9J/T)] (red dashed lines) reached in the 
opposite limit. 



results show that non-magnetic impurities (chain breaks) 
have a quite different effect than paramagnetic impuri- 
ties. In the extreme low-T limit, T <C pv, chain breaks 
will lead to 1/T-behavior with a strength reduced from 
that which would arise from the same density of paramag- 
netic impurities by the fraction of chains of odd length, 
(1 - p)/(2 - p) ~ 1/2. In the opposite limit pv < T 
(but T still <C J), chain breaks lead to a correction to 
the pure result of ss p/[12Tln(2.9J/T)] corresponding 
to an "effective paramagnetic impurity concentration" of 
w p/[3 m(2.9 J/T)] which can be < p for T < J. This 
behavior is illustrated in Fig. 2. 

It has been shown that the Curie-like contribution to 
the susceptibility of Sr 2 Cu03 + 5 can be significantly re- 
duced by annealing under Ar atmosphere indicating that 
this contribution is mostly due to excess oxygen [20] . We 
expect the additional interstitial oxygen ions to be in a 
2 ~ ionization state. Each oxygen ion would then dope 
two holes into the chains. If these holes arc immobile 
they effectively act like chain breaks. It has been pointed 
out that the next-nearest neighbor coupling is not that 
small, J 2 ~ J/16 w 140 K [18]. Whereas we expect 
this irrelevant coupling to cause only small corrections 
to the susceptibility of an isolated chain segment the fact 
that it can bridge a chain break could be significant. A 
first order perturbative calculation, however, shows that 
the correction to the averaged susceptibility is given by 
$Xp ~ Phxp- We can therefore ignore both interactions 
as long as T ^ max(pj 2 , J±) making the model discussed 
here applicable. As shown in the inset of Fig. 3 it is in- 
deed possible to explain the data in [5] for the unannealed 
sample using a 0.6% concentration of chain breaks. If our 
picture is correct the amount of excess oxygen is signifi- 
cantly larger than assumed in [5], 0.06%. 
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FIG. 3: (Color online) Measured susceptibility T(x — Xo) for 
Sr2Cui_ I Pd a; 03+i with x = 0.5, 1, 3% (crosses from bottom 
to top) where a constant xo (see Table I) has been subtracted. 
The green dots represent the best theoretical fit with p as in 
Table I. Subsequent curves are shifted by 3 x 10" 3 . The 
blue solid lines represent (p/4) • (1 — f>)/(2 — p) and the red 
dashed lines Txt,uik + p/[12 ln(2.9J/T)]. Note that fits x = 
Xo+Xijif. +p P ara/(4T) in the regime T > 50 K would be similar 
to the red dashed lines with effective paramagnetic impurity 
concentrations p pala ~ p/[3 ln(2. 9 J/T)] an order of magnitude 
smaller than the nominal Pd concentrations. Inset: Same for 
"as grown" sample of Sr2Cu03+s from Ref. [5]. 



The Sr 2 Cui_ 2; Pd : ].03 + 5 samples used in [7] were not 
annealed. This complicates the analysis of the suscep- 
tibility data because excess oxygen again leads to chain 
breaks in addition to those caused by the non-magnetic 
Pd ions. The amount of excess oxygen hereby seems to 
vary significantly from sample to sample. This is most 
apparent when comparing the measured susceptibilities 
for x = 0.5% and x = 1% which are almost identical 
apart from a small constant shift. We therefore treat the 
impurity concentration p as a fitting parameter. Results 
where we used J = 2200 K and a gyromagnetic factor 
g = 2 are shown in Fig. 3 and Table I. The figure clearly 



TABLE I: Concentration x of Pd ions in experiment com- 
pared to impurity concentration p and constant contribution 
Xo yielding the best theoretical fit. The first line corresponds 
to the "as grown" sample of Sr2Cu03+i from Ref. [5]. 



x (Exp.) 


p (Theory) 


Xo [emu/mol] 


0.0 


0.006 


-7.42 x 10" 5 


0.005 


0.012 


-7.7 x 10" 5 


0.01 


0.014 


-7.5 x 10~ 5 


0.03 


0.024 


-7.5 x 10~ 5 



shows that the susceptibility data are indeed quite dif- 
ferent from what would be expected for paramagnetic 
impurities. Most of the data are taken in the non-trivial 
crossover regime pv ~ T and are well described by our 



theory. Note, that the deviations at low temperatures 
are expected because the suppression of the susceptibility 
due to interchain and next-nearest neighbor interactions 
is not taken into account in our model. For x — 0.005 
and 0.01 the effective impurity concentrations p are larger 
than the nominal Pd concentration, p — x roughly coin- 
cides with the number of chain breaks in the "as grown" 
Sr2Cu03+,5 supporting our picture of additional chain 
breaks due to excess oxygen. For x — 0.03, however, p 
is smaller than the nominal Pd concentration. Possible 
explanations are that for larger doping levels some of the 
Pd ions go in interstitially instead of replacing Cu ions or 
that the sample exhibits some sort of phase separation. 
Future measurements on annealed S^Cui-^Pd^Oa sam- 
ples would be helpful to clarify some of these issues. 

The authors acknowledge valuable discussions with 
W. N. Hardy, K. M. Kojima and G. A. Sawatzky. This 
research was supported by NSERC (J.S., N.L. I. A.), the 
DFG (J.S.) and the CIAR (LA.). Numerical simulations 
have been performed on the Westgrid network. 



[2 
[3 

K 

[5 

[-: 

[8 

[9 

[10 
[11 

[12 

[13. 

[14; 

[15 
[16 

[17 
[18 
[19 

[20] 



J. C. Bonner and M. E. Fisher, Phys. Rev. 135, A640 
(1964). 

S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992). 
S. Wessel and S. Haas, Phys. Rev. B 61, 15262 (2000). 
B. Schmidt, et al., Eur. Phys. J. B 32, 43 (2003). 
N. Motoyama, H. Eisaki, and S. Uchida, Phys. Rev. Lett. 
76, 3212 (1996). 

T. Ami, et al., Phys. Rev. B 51, 5994 (1995). 

K. M. Kojima, et al., Phys. Rev. B 70, 094402 (2004). 

S. Fujimoto and S. Eggert, Phys. Rev. Lett. 92, 037206 

(2004) . 

M. Bortz and J. Sirker, J. Phys. A: Math. Gen. 38, 5957 

(2005) . 

J. Sirker and M. Bortz, J. Stat. Mech. P01007 (2006). 
A. Furusaki and T. Hikihara, Phys. Rev. B 69, 094429 
(2004). 

S. Eggert, I. Affleck, and M. D. P. Horton, Phys. Rev. 
Lett. 89, 047202 (2002). 

I. Affleck, Fields, Strings and Critical Phenomena (Les 

Houches, Session XLIX, Amsterdam, 1988), p. 563, 

edited by E. Brezin and J. Zinn- Justin. 

I. Gradshteyn and I. Ryzhik, Table of Integral, Series, 

and Products (Academic Press, New York, 2000). 

S. Lukyanov, Nucl. Phys. B 522, 533 (1998). 

I. Affleck, D. Gepner, H. J. Schulz, and T. Ziman, J. 

Phys. A 22, 511 (1989). 

S. Qin and I. Affleck, J. Phys. A 32, 7815 (1999). 

H. Rosner, et al., Phys. Rev. B 56, 3402 (1997). 

See EPAPS Document No. [] for the QMC as well as the 

interpolated and extrapolated data. 

From annealed samples we know that the true paramag- 
netic impurity concentration is 0.013% at most [5] and 
can therefore be neglected in the following. 



